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Abstract Analysis of SOHO/LASCO C3 data reveals a discontinuity, inter- 
preted as a shock wave, in plasma density radial profiles in a restricted region 
ahead of each of ten selected coronal mass ejections (CME) along their travel 
directions. In various events, shock wave velocity V w 800-2500 km s _1 . Compar- 
ing the dependence of Alfvcn Mach number Ma on shock wave strength P2/P1, 
measured at R > 10R & from the center of the Sun, to ideal MHD calculations 
suggests that the effective adiabatic index 7, characterizing the processes inside 
the shock front, is largely between 2 and 5/3. This corresponds to the effective 
number of degrees of freedom of motion 2 to 3. A similar dependence, Ma^PiI P\), 
was derived for the Earth's bow shock and interplanetary collisionlcss shock 
waves. All this supports the assumption that the discontinuities in front of CMEs 
are collisionless shock waves. 

Keywords: Coronal Mass Ejections, Initiation and Propagation; Waves, Shock 
1. Introduction 

If the Rankinc-Hugoniot relations are valid at the discontinuity in supersonic 
plasma (or gas) stream parameters, the discontinuity may be interpreted as a 
shock wave. Even to roughly estimate the validity of these relations for a shock 
discontinuity driven in front of fast CMEs is no easy task. This is first of all 
associated with the fact that a shock front is difficult to identify and register 
in the coronal conditions. Nevertheless, several such attempts for individual 
events have been made in IVourlidas et all 120031 [Eselevich and Es clevichl I2007| 
IManchester IV et aU [2008: Ese levich and Eselevichl[2H051IOntiveros and Vourlidasl l 
120051 Thereafter Paper 1). 

Two different approaches may be singled out which have been applied to 
identifying a shock front on coronal images. The first one is numerical simulation 
within the framework of the ideal magnctohydrodynamics (MHD) of the entire 
process of CME formation and propagation. In this case, the manner in which the 
CME is formed may be arbitrary. The results of such a simulation were compared 
to the experiment in order to identify the shock front ( |Vourlidas et aL,~2 003; 
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[Manchester IV et al., 2008] ). Note that in this case one can directly estimate the 
relation between thermodynamic parameters at the discontinuity because they 
are included in the model. 

The second approach relies on the fact that there is no shock discontinuity 
if the velocity of the generating CME remains below the critical velocity, uc', 
while such a discontinuity does exist at supercritical velocities. Regardless of 
the difficulty in estimating the critical velocity in the corona, this approach 
enabled us to demonstrate that the transition does exist, and consequently the 
observed discontinuity is a shock wave ( Eselevich and Eselevich, 2008[ ). There is 
some difficulty in identifying the shock fronts in coronal images stemming from 
the fact that the CME frontal structure and the shock front in front of it are 
identical geometrically and are often rather close to each other. We employed two 
different complementary methods to make sure that the observed discontinuity 
in the brightness distribution is the shock front rather than the CME frontal 
structure behind it: 



1. We examined the conditions in the disturbed zone in front of CMEs with 
different velocities u. For CMEs with velocities exceeding u, the shock- wave- 
related discontinuity was observed in the frontal part of a rather extended 
disturbed zone. 

2. The processes of the disturbed zone evolution and shock wave formation were 
investigated for some CMEs, in the coordinate system associated with their 
frontal structures, when the CME velocity u passes the critical velocity uc 
( [Eselevich, 20100 . 

Note that the value of the critical velocity at which a shock wave formed ap- 
peared to be rather close to existing estimates of the Alfven velocity in the 
corona ( |Mann et al., 19 99). These methods allowed the shock front width to be 
measured for several tens of CMEs at various distances (up to 6Rq) from the 
center of the Sun (i? Q is the Sun's radius). 

Thus, whether by one or another method, we can identify a shock wave in 
the corona; so the next obvious step would be to try to verify the validity of the 
Rankinc-Hugoniot relations at the shock. We are restricted, however, in that we 
cannot measure just any plasma parameter in the corona. White-light coronal 
images allow the shock discontinuity velocity to be determined (from the change 
in the discontinuity position with time) and the density ratio estimated at the 
discontinuity, relying on some assumptions. Apparently, it is these values that 
should be taken as a basis for the analysis. Moreover, similarly to the Earth's 
bow shock, we may apply simple relations of the ideal MHD at a flat shock 
discontinuity to investigate shocks in the corona, i.e. make use of the assumption 
that the shock front width is smaller than the radius of its curvature. 

The aim of this work is to obtain Rankinc-Hugoniot relations at the flat 
shock discontinuity, within the framework of the ideal MHD, and to estimate 
their validity for a number of identified shocks driven ahead of CMEs. 
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2. Data and method of analysis 

The analysis involved LASCO C3 coronal images with a 4 to 30 Rq field of view 
(Brucck ner et al., 1995[ ). The images contain many stars that are much brighter 
than the corona. Such point objects as stars can be easily removed by a sigma 
filter. The filter is applied to calculating the mean value and rms deviation in a 
square neighborhood centered at each image pixel, excluding the pixel itself. If 
the brightness of the central pixel exceeds the mean value by more than a certain 
value, it is replaced by the mean. The filtering of all images involved a 15x15 
pixel neighborhood and a maximum deviation equaling three rms deviations, 
while the filtering procedure was iterated ten times. After filtration, the images 
of bright stars are effectively eliminated, the coronal brightness signal remaining 
almost unchanged. 

The images thus processed were presented as difference brightness AP = 
P(t) — P(to), where P(t) is the coronal brightness at instant t, corresponding to 
the event in question, P(to) is the undisturbed brightness at fixed instant to cho- 
sen well before the event. Thus, the background brightness due to the stationary 
distribution of plasma in the corona was eliminated. Calibrated LASCO images 
with full brightness P(t) (in units of mean solar brightness, Pmsb) enable us to 
use the difference brightness to estimate the mean density of the matter causing 
the change in brightness, if we make an assumption with regard to its spatial 
size along the line of sight. 

The difference brightness was used to examine the dynamics of CMEs, as well 
as of the disturbed zone and shock wave in front of them. It has proved to be 
more convenient to employ maps of difference brightness isolines rather than 
the traditional grey-scale images for selecting and generally analyzing events. 
The arrangement of isolines provides for a better understanding of the character 
of the difference brightness distribution, direction and value of the gradients. 
Isolines levels were selected depending on the event and the distance at which it 
was considered. 

It is more convenient to make quantitative measurements using difference 
brightness profiles drawn in a certain direction. For this purpose, we used radial 
profiles AP(R) of difference brightness drawn from the center of the Sun along 
the direction of the fixed position angle PA as well as nonradial profiles AP(r) 
constructed from a specified point and in a specified direction. The notations 
differ in distance R, plotted from the center of the Sun, and in distance r, plotted 
from another point on the image. All the distances in the plots and figures are 
in units of one solar radius (Rq). In some cases, difference brightness profiles 
correspond to the finite angle average. The averaging had a purpose of improving 
the desired signal-to-noise ratio. 

3. Relations at a flat shock discontinuity in the ideal MHD 
appr oximat ion 

Relations at a discontinuity link various plasma parameters on either side of the 
shock discontinuity. They enable us to relate the Mach number in the medium 
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where the discontinuity moves to the parameter ratio at the discontinuity. Under 
the coronal conditions (at least, at l-3Oi?0 from the center of the Sun), the Alfven 
velocity va far exceeds the sound velocity cs- We can therefore take va for 
the typical velocity of disturbance propagation and find the dependence for the 
Alfven Mach number Ma = u/va, where u is shock wave velocity. An expression 
for Ma at the flat shock discontinuity may be obtained in its analytical form 
(its derivation is presented in Appendix A). 

It is so far impossible to measure magnetic field and temperature jumps at the 
shock discontinuity in the corona. But we can estimate the density ratio pi / p\ at 
the discontinuity from the white- light coronal brightness (index 1 means that the 
value corresponds to the undisturbed region ahead of the shock; 2 to the region 
behind it). Thus, our prime interest will be with the Ma^PiI pi) dependence. The 
angle between the magnetic field direction ahead of the front and the normal to 
the wave front 0Bn, adiabatic index 7, and the gas to magnetic pressure ratio 
/3i still remain free parameters in the equation. 

Figure 1 presents sets of Ma[.PiI P\) curves for 7 = 2, 5/3 and 1.46 at different 
0Bn and (3\. The left panel in Figure 1 shows the Ma(p2/ Pi) dependences for 
(3i = and Osn = 90° and 45° (solid and dashed lines respectively). In the right 
panel Figure 1,0 Bn = 90° and /3i = and 1. The dependences Ma(p2/ 1 Pi) at 7 = 
5/3 and /3i = coincide with similar dependences in IKantrowitz and Petschecl 
119661 as presented in IKennel et al\ 119841 Figure 1 implies that the form of 
Ma^PiI P\) for quasi-perpendicular shock waves (0Bn > 45°) depends little on 
Obu and (3\ (at < /?i < 1), but varies considerably depending on 7. 



4. Analysis of experimental data 
4.1. Event selection criteria 

The main criterion for selecting events for the analysis was a reliable identifica- 
tion of shock fronts in the events. The identification was based on the following, 
previously obtained, results ( |Eselevich, 2010| ): 
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1. Ahead of CMEs with velocities exceeding w 700 km s -1 , it is possible to 
record a shock wave and correctly measure its front width Sf based on Mark 
4 and LASCO C2 data. Our measurements have also shown that the shock 
front width 8p is of order of the free proton path at 1.2i? Q < R < 6Rq, at 
least in a limited region located in the direction of the CME propagation. 
This provides evidence for the collisional mechanism of energy dissipation in 
the shock front at these distances. 

2. A new discontinuity with width S F -C Sf is observed to form at R > 10-R© in 
the frontal part of the shock front. Within the experimental error, the value of 
S F w O.l-O.2i? is independent of R, being defined by the LASCO C3 spatial 
resolution. This brightness profile transformation from the front of width Sf 
to the discontinuity of width 5 F -C Sp was interpreted as transition from a 
collisional to a collisionless shock whose front width could not be resolved in 
coronal images. 

The authors of Paper 1 registered brightness distribution discontinuities inter- 
preted as shock waves, at the front of 13 selected CMEs moving at over 1500 km 
s _1 . Since the heliocentric distances at which the discontinuities were observed 
generally exceeded 6-10i?Q, and the minimum width of these discontinuities 
did not exceed 0.1-0.2i? Q , they appear to be associated with collisionless shock 
waves. 

The flat shock wave condition S f <^Rf(Rf^s the radius of wavefront curva- 
ture) is better satisfied at large distances because it is there that fronts of width 

0. 1.0.2i? Q are encountered. Moreover, the shock wave strength is sufficiently 
great {pi/ ' pi > 2) at these distances; i.e. the shock wave has already been formed 
and is stable. Therefore, we chose for our analysis those instances of time when 
the shock front was farther than lOi?0. 

Comparing the experimental data with the calculations of relations in the 
shock front requires absolute measurements of density; it is therefore very im- 
portant to select events with minimised influence of factors affecting the mea- 
surements. We have chosen ten CMEs having shock waves in front, while keeping 
in mind the following requirements: 

1. The origin of the CMEs was near the limb and they propagated in the plane 
of the sky; i.e., their measured positions and velocities were close to the true 
values. 

2. 12-24 hours before each of the selected events, there were no other CMEs 
capable of noticeably changing the undisturbed solar wind conditions in the 
coronal region under study. 

These events are listed in Table 1. The third column in the Table lists the 
coordinates of the CME-associatcd source region on the disk. The coordinates 
for the first two events were borrowed from ICremades and Bothmeri 120041 for 
the others from the catalogue lYashiro et ail 120041 Since the origin of all the 
ten CMEs was near the limb, let us call them 'limb CMEs'. According to the 
classification accepted in the catalogue lYashiro et al.i 120041 they include halo 
(events 3, 4, 7, 8 in Table 1) and partial halo CMEs (event 2). Let us not 
confuse limb CMEs with coronal ejections appearing at longitudes below 60-70°, 
which arc therefore no limb CMEs. 
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Figure 2. An example of a CME moving at a supersonic velocity (difference brightness 
isolincs, LASCO C2 data). This figure includes the notations employed for referencing and 
constructing difference brightness profiles required for measurements in the shock front. The 
relevant explanations are given in the text. 



4.2. Registering shock fronts ahead of CMEs 

The shape of the CME outer boundary (and the adjacent shock wave) is usually 
nearly circular when the CME velocity exceeds the critical value uq (Figure 2). 
The center of this circle (point in the figure) may be defined by distance 
Rc from the center of the Sun (point S) and on the position angle PAc- The 
difference brightness profiles AP(r) used to determine shock wave position were 
drawn from the center of the circle. The position of each profile was given by angle 
a to the direction of CME propagation. This angle is positive counterclockwise. 

A shock associated discontinuity in difference brightness distributions is gener- 
ally observed in a limited region in the direction of CME propagation (Esclcvi ch, 2010[ ).[ 
The value of a at which shocks were recorded did not exceed ±10° in all the 
ten selected events. For this range of angles a, the distances from the center 
of the Sun, R, and from the CME center, r, to the shock front may be related 
by a simple relation R w Rc + r, with an accuracy of O.li?0. However, even 
a-few-degree change in a and a-few-pixel change in the CME center position 
significantly affect the discontinuity width as recorded in difference brightness 
profiles. An accurate choice of the O center position almost always provides the 
minimum discontinuity width comparable to the LASCO C3 spatial resolution. 

Difference brightness profiles AP(r) drawn from the CME center were used 
to register shock fronts and measure their positions and amplitudes. Figures 3-6 
present examples of such distributions for four events in Table 1. Each of the 
figures corresponds to one of the events, and each plot in the figures corresponds 
to a certain instant of time. The time difference used to construct the difference 
brightness, the CME center position (angle PAc and distance Rc), and the 
direction along which the profile was drawn (angle a) are given for each plot. In 
all the cases, we used averaging in the range of angles da — 2°. For convenience, 
each of the plots is reconstructed depending on the distance R from the center 
of the Sun, though the scanning direction may slightly differ from the radial 
direction. Each plot presents two distributions: 1) for the instants when CMEs 
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and shocks were recorded and 2) the difference brightness distributions for two 
neighboring instants before the CMEs appeared in the C3 field of view (the 
distributions are denoted by black and light circles respectively). The latter 
make it possible to estimate the noise level; in particular, they were used to find 
the error (lcr) in measuring the shock wave amplitude 5P. At large distances, 
the noise level is comparable with the amplitude in some of the events. 

All the events have common features, but differ in some respects. The dif- 
ferences require some commentary. The upper plots in Figures 3 and 4 show 
the stage when a collisional shock wave with front width Sf transforms into 
a collisionless discontinuity with front width 8* F -C Sf, for the 1999 July 25 
and 2007 December 31 events, respectively. Such a transformation usually takes 
place at R < IORq. The location and amplitude of 5P were determined only for 
shock discontinuities with front width S F at > 10i? Q . Another event occurred 
in the same coronal region approximately 12 hours before the CME, during the 
1999 July 25 event (Figure 3). As a result, the undisturbed difference brightness 
constructed for a slightly earlier instant of time (light circles in the two upper 
plots) has a mean value that is significantly different from zero. In determining 
5P, we used averaging over a region behind the front whenever possible (Figures 
3-6). The front width 8* F slightly varies but is almost constant at >10-15i?o (the 
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minimum value is 0.1-0.2i?Q), and the observed variations may have been due 
to noise. 

The plots in Figures 3-6 as well as similar plots for other events were used 
to determine the shock wave distance from the center of the Sun (Column 4 in 
Table 1) and shock wave amplitude (Column 5) at a given instant of time. The 
obtained distances were employed to calculate shock wave velocities (Column 6). 
The velocity calculations took into account information on shock wave locations 
at < lOi?© at earlier instants of time, not included in Table 1. The velocity at a 
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given instant of time is the average velocity for the intervals before and after this 
instant, except for the last instant, for which the velocity is simply the velocity 
as determined for the previous interval. 

4.3. Estimating shock wave strength and comparing the results with MHD 
calculations 

4-. 3.1. Experimental values for estimating the shock wave strength 

Our aim is to compare the experimental dependence Ma{p2/ Pi) with the anal- 
ogous dependence derived from MHD calculations. The shock wave strength is 
defined as P2/P1 ~ l + SN/N , where SN is the absolute jump in electron density 
in the shock front, Nq is the undisturbed plasma electron density immediately 
ahead of the shock front. To find Ma, we have to know the local Alfven velocity 
(Va) and the velocity (Vsw) of the solar wind relative to which the shock waves 
move. 

Since iVo, Va and Vsw cannot be determined directly from the available 
experimental data, let us use their mean values typical of the quasi-stationary 
slow solar wind (SW) flowing in the streamer belt. The distance R dependences of 
Vsw, Va, and iV (Figure 7) were taken respectively from QWang et at, 2000] ), 
QMann et al, 1999[ ), and flSaito et al, 1977[ ). The validity of using the above 
dependences is obvious enough because the shock front parts under analysis 
are located near the axis of the CMEs propagating in the slow SW. This is 
also confirmed by the successful application of the dependences Vsw(R) and 
Va(R) in the analysis in (Esclcvich and Esclcvich, 2008; sec Figure 3) and, to 
some extent, of the dependence N {R) in Paper 1. Applying No(R), however, 
has a limitation mentioned above - there should be no other CMEs capable 
of changing the undisturbed SW conditions considerably in the region under 
study in the previous 24 hours. Column 7 in Table 1 gives calculated values 
of Ma — {V — Vsw)/Va, where Vsw and Va values derived from the above 
experimental dependences. 

The density jump SN in the shock front remains an unknown value and 
is to be determined. It is proportional to brightness jump 8 in the difference 
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Figure 7. The experimental distributions employed to calculate Ma and P2I P\- The upper 
panel: the solid curve is the Alfven velocity Va{R) in lMann et aU J 19991 the dashed curve is 
the slow SW velocity Vsw(R) m Wang et al. 2000- the lower panel is the electron density in 
the corona in ISaito et qi.llT5771 



brightness distributions (Figures 2-6). The relationship between S and 6N de- 
pends on the unknown plasma distribution in the front. However, by specifying 
mean line-of-sight size I of the shock front, it is possible to calculate the mean 
density jump from the measured brightness discontinuity, using known relations 
JBilling, 1966] ) . As long as we are dealing here with limb events, we may assume 
that the shock front is in the plane of the sky. 

In Paper 1, the mean size of the shock front - I = IRq - was used to calculate 
SN from 5P for all the cases of the nonlimb CMEs. In what follows we will try 
to justify quantitatively a somewhat different value of I for the limb CMEs in 
Table 1, common for all the events under investigation as well. 

4-3.2. Estimating the ultimate compression in a collisionless shock wave front 

In the front of a quasi-perpendicular (Obu > 45°) collisionless shock wave, heat- 
ing of ions can occur in a plane perpendicular to the magnetic field direction 



(Zimbardo, 2009). This transverse motion is two-dimensional, which corresponds 
to two degrees of freedom (i — 2) and consequently to the adiabatic index 
7 = (i + 2)/i = 2. In the case of isotropic turbulence in the shock front, the 
number of degrees of freedom will tend to i = 3, and the adiabatic index to 
7 = 5/3 (Sa gdeev, 1966| ) . It is therefore reasonable to expect that the adiabatic 
index 7 will be between 5/3 and 2 in the front of a quasi-perpendicular colli- 
sionless shock. This means that maximum pij p\ values at large Mach numbers 
must not exceed 4 (see Figure 1 or Appendix). 

This assumption is supported by the results of direct measurements in the 
Earth's bow shock and interplanetary shocks. These measurements do not pro- 
duce pij P\ > 4 even for the largest Mach numbers. To confirm this, the up- 
per panel of Figure 8 presents the experimental dependences Ma{p2/ 1 Pi) for 
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Interplanetary shack 

Irom Zastenkei el al., 
1 -29.01.1981 
2-05.03.1981 

3 - 12.03.1981 

4 - 18.04.1981 

5 - 10.05.1981 

6 - 15.05.1981 
7-25.07.1981 

Earth's bow shock ■ 

from Formisano el al., 
1 - 12.02.1969 



9- 06.11.1977 

10- 27.11.1977 
11 - 12.12.1977 

12- 13.12.1977 

Irom Bale era/., 2005 

13- 25.12.2000 



CME's shocks 

■ - 20.04.1998 
A - 20.09.1997 
V - 25.07.1999 

♦ - 21 .04.2002 
□ - 16.08.2002 

• 04.11.2003 
l> - 23.08.2005 
► - 05.09.2005 
+ -31.12.2007 
A - 25.03.2008 



Figure 8. The experimental dependence Ma(P2/ pi)'- the upper panel is the bow shock (dark 
squares) and interplanetary shocks near the Earth's orbit (light squares); the lower panel 
displays ten shocks associated with limb CMEs. The solid and dashed curves denote ideal 
MHD calculations for perpendicular (Osn — 90°) shocks at 7 = 5/3, 2 and 3. 



bow (solid marks) and interplanetary (light marks) shocks based on data of 
various authors. The values of pil p\ in these dependences were taken from 
( Zastenker et al, 1983 ) and ( Formisano et al, 197 3; Greenstadt et al, 1980 Bam e et al, 1979[ | 



Bale et al, 20051 for interplanetary and bow shocks respectively. In some cases, 



the calculation of the Alfven Mach number was based on hour-averaged OMNI 
data (http://omniweb.gsfc.nasa.gov/) on the magnetic field, SW density and 
velocity (for the bow shock). All the obtained p^j 'pi do not exceed 4. This is 
also consistent with the Mms (P2/P1) plot in Figure 6 in lFormisano et 00119731 
That figure implies P2/P1 < 4 for all 42 crossings of the bow shock with large 
magnetosonic Mach numbers Mms — 4 — 12 (by definition, Ma > Mms)- 

The upper panel of Figure 8 implies that the experimental points are mainly 
located near the calculated curves for 7 = 5/3 and 7 — 2, i.e. in the collisionless 
shock front 7 > 5/3. Note that the effective adiabatic index 7 is smaller in the 
undisturbed SW. For example, ITotten et al\ 119951 derived the empirical, mean 
adiabatic index for protons - 7 = 1.46 =b 0.04 < 5/3 - using Helios 1 data for the 
SW moving at V S w = 300 - 800 km s" 1 . The fact that shocks arc characterized 
by greater values oi7«5/3 — 2 must be a result of collisionless processes inside 
the shock front. 
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Table 2. 



a, deg. 


R, Rq 


c t~) wir\ — 13 n 

6P, xlO lo Pmsb 


-20 


24.8 


2.78 


-17 


26.0 


2.86 


-14 


26.7 


3.02 


-11 


26.9 


2.00 


-1 


25.4 


3.24 





25.5 


3.11 


+3 


26.4 


2.57 


+6 


25.7 


3.37 


+10 


26.8 


2.94 


+18 


26.4 


3.36 


+19 


26.3 


2.65 




26.1 + 0.66 


2.9 + 0.4 



4--3.S. Determining the effective line-of-sight size I of the shock front 

By analogy with the heliosphere it is reasonable to expect that the adiabatic 
index 7 > 5/3, and maximum density ratio is less than 4 for a collisionless shock 
in the solar corona. Let this conclusion serve as a basis for estimating I. For the 
purpose, we will examine the 2002 April 21 CME in more detail (number 4 in 
Table 1). This CME had the highest velocity of the ten selected events, and a 
maximum compression was observed in the shock front in front of the CME at 
about 26i?0 from the center of the Sun at £=03:16:57. To enhance the accuracy, 
we constructed difference brightness distributions for several angles a in the 
range —20° to +20°. Five of the distributions are shown in Figure 9. Besides, 
to determine the brightness jump SP the brightness profile was averaged on 
either side of the discontinuity. The event was accompanied by an intense flux 
of energetic particles which manifested themselves as bright dots and scratches 
on the coronal images. In some cases, filtration failed to completely remove 
such noise from the images. However, it can easily be detected in the brightness 
profiles (one of the examples is given in the upper plot of Figure 9). These parts 
of the brightness profiles were not used to determine SP. 

Table 2 summarizes the measured 5P values for various angles a and the 
distances at which they were observed. The distances differ slightly, indicating 
that the shock front is inhomogeneous along the direction of angle a. Deter- 
mination of the front velocity and density discontinuity employed their mean 
values. The mean brightness discontinuity was SP = 2.9 x 10~ 13 Pmsb, and the 
mean velocity V = 2422 km s _1 . This velocity corresponded to Ma = 7.94 
at the average distance R = 26. IRq, where the shock front was recorded. The 
undisturbed density at the distance was No = 1.3 x 10 3 cm~ 3 . 

The brightness jump SP at a given distance R allows us to calculate SN and, 
taking Nq into account, derive pij P\ from it. As was stated above, we must 
set the shock front size I along the line of sight for the purpose. The curve 
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PA C =272°, R C =14.1R 

• -03:16:57 -00:16:58 
0-00:40:50-00:16:58 

a = +18° 




Figure 9. The difference brightness AP(R) distributions for several angles a for the CME 
of 21 April 2002; the shock front ahead of this CME had the highest compression of all the 
events under study. 



in Figure 10 demonstrates how the calculated pij P\ varies with the selected 
size I (in units of solar radius) under the conditions of the event under study 
(SP = 2.9 x 1(T 13 P A /,s,b, R = 26. 1R Q , N = 1.3 x 10 3 cm- 3 ). The main critical 
parameter is the adiabatic index 7 because the /3 and Obu dependences are rather 
weak. The range of maximum pij p\ (i.e. as Ma — > 00) for 7 between 5/3 and 2 
is indicated by hatching in Figure 10. 

Let us assume that an extreme case is observed in the event: 7 = 5/3, f3 = 
and Qbu = 90°, i.e. there is a perpendicular shock with maximum compression. 
Then, the MHD calculations imply that the shock wave strength pij P\ = 3.74, 
being close to the limiting value of 4. The plot in Figure 10 indicates that this 
strength corresponds to I = 6.5Rq. This size was employed to calculate 8N and 
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Figure 10. The pil p\ values depending on the effective length I along the line of sight. They 
were calculated for the shock wave parameters at 03:16:57 ahead of the CME of 21 April 2002. 

Pil P\ m a ll the ten selected events, their values are listed in columns 8 and 9 of 
Table 1. 

4-.S.4- Experimental dependence Ma(P2/ ' Pi) for shocks in the corona 

The calculated values in the lower panel of Figure 8 are used to compare the 
experimental dependence Ma{j>-xI P\) for ten shocks with calculation results for 
perpendicular shocks (Obu = 90°) at 7 = 5/3 and 7 = 2 (solid lines). The error 
bar due to the error in determining the brightness jump is indicated for the 
experimental points. 

The dependences Ma(p2/ Pi) suggest that as a whole they fit the calculated 
curves reasonably well in spite of the considerable scatter of experimental points. 
This allows the following conclusions to be made: 

1. The choice of the mean value I = 6.5Rq has proved to be reasonable enough 
despite the fact that the shock fronts in question were observed at various 
distances (10R & < R < 3Oi? ) and had velocities V « 800-2500 km s _1 . 

2. The comparison results agree with the assumption that the effective adiabatic 
index 7. characterizing the processes within the front, is largely between 2 
and 5/3, which corresponds to the effective number of degrees of freedom of 
motion being 2 to 3. This also testifies that the structures under investigation 
are collisionlcss shock fronts. 

Of special notice are several outstanding points which correspond to shocks 
with sufficiently large Mach numbers but with small pil P\ (outlined in upper 
and lower panels of Figure 8). They all lie to the left of the calculated curve with 
7 = 3. This corresponds to the effective number of degrees of freedom \ 1. If these 
points do not result from errors, we have to admit that in these cases the shock 
is not in a stationary state (physically justified dependence Ma(P2/ 1 Pi) becomes 
invalid for the points). It may result from various reasons. For example, in the 
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Figure 11. (light circles) and P2/P1 variations (dark circles) with distance in the events on 
31 December 2007, 25 July 1999, and 16 August 2002. 



case of the bow shock, the nonstationarity of the front may result from a change 
in the parameters of the upflow undisturbed SW. 

The experimental results in Table 1 allow certain conclusions to be made with 
regard to how the shock front propagates in a coronal plasma where concentra- 
tion rapidly decreases with distance. On the example of three different events in 
Figure 11, one can see that the shock wave strength and Alfven Mach number, on 
the average, either change insignificantly or increase with distance at distances 
up to R < 3Oi?0. Such behaviour is also typical of other events. 



5. Comparison with results of Paper 1 

Density ratios pij p\ (marked Tqr) have been found in Paper 1, for eleven shocks 
driven by CMEs moving at 1500-2000 km s _1 . The parameters of the events are 
in Tables 1 and 2 in Paper 1. They are numbered 1, 2, 4, 6, 7, 8, 9, 10, 11, 12, 
14. The table velocity and distance values were used to calculate Ma values. 
The asterisks in Figure 12 mark the values corresponding to the events on the 
Ma{p2/ Pi) plot. Since the pij p\) values were derived in Paper 1 for I = IRq, 
the density ratios were recalculated for I = 6.5i?0 to compare them with our 
findings. The recalculated values are marked by dark circles in Figure 12. The 
light circles indicate our values (the same as in the lower panel of Figure 8) . 

It is obvious that there is a significant difference between findings in this paper 
and Paper 1: the density ratios in Paper 1 are much smaller than the values we 
obtain here. Note that the method for determining pij p\ is identical in both 
papers. The only difference is the value of IRq in Paper 1 and 6.5-Rq in our 
paper. There may be several reasons for the difference. We list them in the order 
of their influence on the result : 
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Figure 12. The experimental dependence Ma{P2/ Pi) constructed from the P2/P1) values 
derived in this paper (light circles) and in Paper 1 for I = IRq (asterisks) and recalculated for 
I = 6.5i?0 (dark circles). 

1 . presence of other CMEs capable of markedly changing the undisturbed solar 
wind density in the coronal region under study, less than twelve hours before 
the CME (events 1, 2, 4, 8, 10, 12 in Paper 1); 

2. a non-stationary shock front, when the registration point is at R <1O-I5i?0, 
and p 2 /pi < 2.0 - 2.5 (events 1, 4, 6, 8, 9, 10, 11, 12, 14); 

3. the CME source is rather far from the limb (nonlimb event) resulting in a 
noticeable error in determining distance R at which the shock is registered 
(and thus p\ and Va) and its velocity V . This may also result in decreased 
intensity of scattered light and, thus, decreased measured value of S (events 
1, 2, 6, 8, 9, 10, 11, 12). 

Finally, the last comment. Paper 1 does not include pil Pi > 2.8. But if we set 
/ = IRq for our events with pil Pi > 2.5 at I = 6.5Rq, we will obtain P2I Pi > 
10. Obviously, this value far exceeds the maximum values of the density ratios 
observed for collisionless shocks in the heliosphcre (upper panel in Figure 9). 

6. Conclusions 

1 . A brightness discontinuity, interpreted as a shock wave, was demonstrated to 
be registered in a limited region in front of each of the ten selected CMEs 
(propagating in the plane of the sky) along their travel directions. In various 
events, the shock wave velocity was V « 800-2500 km s _1 . 

2. Comparing the dependence of the Alfven Mach number Ma on the shock 
wave strength pi I Pi, measured at R > 10i?© from the center of the Sun, 
with ideal MHD calculations suggests that the effective adiabatic index 7, 
characterizing the processes inside the front is largely between 2 and 5/3. 
This corresponds to the effective number of degrees of freedom of motion 
2 to 3. A similar dependence Ma(p2/ Pi) was obtained for Earth's bow and 
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interplanetary collisionless shock waves. All these substantiate the assumption 
that the discontinuities in question in front of the CMEs at R > IQRq are 
collisionless shock waves. 



Appendix 

A. Deriving the dependence of the Alfven Mach number on the 
density ratio at the fiat shock using the ideal MHD 
appr oximat ion 

A.l. Notations and initial equations 

Let us consider a flat shock front. Index 1 will denote the undisturbed plasma 
region ahead of the front, index 2 the region behind the front. A front-associated 
coordinate system is used, i.e. the plasma velocity ahead of the front is the same 
as the front velocity. The front surface is in the yz plane, and the x axis is 
perpendicular to the surface. The coordinate system is chosen such that the 
magnetic field vectors and velocities lie in the xy plane. The magnetic field 
vector B\ is at an angle, 0Bm to the x axis (i.e. to the normal to the front). 
Thus, the case of 9bu = 90° corresponds to the perpendicular shock when the 
magnetic field vector is in the plane of the front. 

According to [Priest! 119821 the following conditions are satisfied at the shock 
discontinuity (in the CGS system): 

P2V2x = PlVlx (1) 

p 2 + B|/8tt - BlJA-K + p 2 v\ x = V x + B 2 /8n - BfjAn + Pl v 2 lx (2) 

P2V2xV2y ~ B 2x B 2y / AlT = plV lx Vi y - B lx B ly / AlT (3) 

( P2 + BI/8tt)v 2x - B 2x (B 2 • v 2 )/4tt + (p 2 e 2 + p 2 v 2 2 /2 + B\l%v)v 2x = 

(pi + Bj/8ir)v lx - B lx (B 1 ■ vi)/4tt + (pid + pivf/2 + Bf/8ir)vi x (4) 

B 2x = B\ x (5) 

v 2x B 2y — v 2y B 2x = v\ x B\ y — v ly B lx (6) 

here v is velocity, B is magnetic field, p is pressure, p is density, e is internal 

p 

energy defined by e = — for the ideal polytropic gas, where 7 is the 

adiabatic index. 

Equation (1) is the mass conservation condition; equations (2) and (3) are the 
conservation conditions respectively of the x and y components of the momen- 
tum; equation (4) is the energy conservation condition. Equation (5) (following 
from div B = 0) is the conservation condition of the normal component B x of 
the magnetic field. Equation (6) is the continuity condition of the [v x B] value 
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following from the electric field tangential component continuity and magnetic 
field frcczing-in condition. 

Equations arc written here in the general form, but the velocity and magnetic 
field z-componcnts arc zero throughout the selected coordinate system. At the 
shock discontinuity, the velocity changes - only the x component of velocity is 
present in region 1, whereas the y component is present in region 2 as well. The 
magnetic field, namely its y component, also varies in the front, whereas the a; 
component remains unchanged. 



A. 2. Derivation and solution of main equations 



Let us introduce the notation for the ratio of tangential (about a wave front) 
magnetic field components: 

h = B 2 y/Biy, 

and the notation: 

r = 7/(7-1). 

We will additionally use the following expressions: Va = B/^Airp for the Alfven 

p 

velocity, Ma = v\ x /va 1 for the Alfven Mach number, and j5 = for the 

B /o7T 

gas to magnetic pressure ratio. 

Equation (3), with due account for (1) and (5), may be used to express the 
tangential velocity component appearing behind the front: 

v 2y = cos 2 6 B n tan 6 Bn 

4npiv lx 

Substituting the v 2y expression into (6) and taking (1) into account, we get: 
v 2x Pi 1 f h-1 2 a \ rvs 

= — = T (1 + -TT2- C0S °Bn) (7 

Vi x p 2 h M\ 

The expression (7) relates the density ratio p 2 /pi and the Alfven Mach num- 
ber Ma- However, this expression also includes the ratio between magnetic field 
tangential components, h; thus we must construct and solve an equation system 
to express Ma in terms of h. 

Dividing both sides of Equation (2) by p\v\ x with due account for (5) and 
expression for v 2x /v\ x yields: 

1/ h — 1 o , , , Si (f - h 2 ) sin 2 9 Bn 

- h (1 + cos 2 ^) (1 + a 2 ) = ^ + + 1 (8) 

We introduce the notation a 2 = 5- here. 

P2V 2x 



SOLA: hugo-re.tex; 13 April 2010; 0:58; p. 19 



Esclcvich and Esclcvich 



3 

Dividing both sides of Equation (4) by — ^-^ and substituting expressions for 
V2 y and V2x/vix, we obtain a second equation: 

, 1 , h-l 9 „ n9 2{h- 1) sin 2 B Bn 
^ + 2Fa2) 5? (1 + ~M\ C ° S W + Ml 

(h-\f^2Q Bn _rft 

Thus there are two equations to express Ma in terms of /i. Expressing ct2 
from the first equation and substituting it into the second one produces, after 
some simplification: 

1 [{m 2 a - cos 2 e Bn ) (m|(i + /i - 2r) + fy? x r -{h- i)(2r - 1) cos 2 e Bn ) 
+h sin 2 e Bn (Af^(/i(r - 2) + r) + (h - i)0(r - 1) + r) cos 2 6> Bn )) = o 



hM 



This equation has two solutions for M\. 

Ml = 2{l + l_ 2T) (-^i r + 2(1 + (h - 2)r) cos 2 e Bn 

-h(h(T - 2) + T) sin 2 6 Bn 
±^Jh?{^T - l) 2 cos 4 9 Bn + (ftT + (/i(r - 2) + T) sin 2 e Bn y 

- 4(r - 1) cos 2 Bn (/3ir + (i + h 2 - (i + /i)r) sin 2 e B „) ) ) 



(10) 



(11) 



The solution with '+' before the root gives the trivial value M\ = at 
0Bn = 90° and physically infeasible values Ml < at 9 Bn ^ 90°. Of interest is 
therefore the second solution with ' — ' before the root. 



A. 3. Extreme cases 



Let us consider some extreme cases of the solution. We may set /3i = for plasma 
totally controlled by the magnetic field, and 7 = 5/3 and correspondingly V = 2.5 
for gas with three degrees of freedom. 

In the perpendicular shock case (6 Bn = 90°), the expression for Ml is 
simplified and takes the form: 

2 _ ft(2-5 + 0-5/1) 

' A 4-h 

If Ml = 1, then h = 1, and if Ml — > oo, then h — > 4. As is evident from (8), 
the density ratio in the shock front in this case simply equals the magnetic field 
tangential component ratio, i.e. pij P\ = h. 
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In the parallel shock case (0g n = 0°), M\ = 1 for any h because it is an 
extreme case of shock wave unaffected by the magnetic field. 

In the general case of f3\ ^ 0, M\ > 1 for h = 1 because the characteristic 
velocity of disturbance propagation in the plasma will be the fast magnetosonic 
wave velocity rather than the Alfven velocity va ( [Priest, 1982[ ) 



where c s is the sound velocity. In this case it is more correctly to try to find the 
dependence of the magnetosonic Mach number Mms = v/vms on density and 
magnetic field ratios. 
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